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Abstract 

Tne  dispersion  curves  lor  surface  plasma  wr.ves  (SPW)  in  An  have  been  dotern.ineu  iron,  calcu- 
lated reflectivity  minima  as  exhibited  by  attenuated  total  reflection  (ATR  spectra)  lor  the 
p.-:  - ;.»-a  i • l ai  cor. J'igura tio:.  and  from  the  direct  calculation  of  tlw  dispersion  rc’arion 

rcr  the  same  configuration.  The  dispersion  curves  for  Au  have  been  determined  by  meusurii.ft  t .*•« 

ATrt  opt ctr*»  lor  the  prion-dielectric-metal  (PD.v»)  configuration,  by  calculating  the  ATP.  spectra 
frcir.  puoliBhC'l  optical  constants  and  from  the  direct  calculation  of  tne  dispersion  relatic:  for 
tne  Pis),  configuration.  Vc  have  found  two  general  types  of  solutions  from  the  direct  calc  latior. 
ox  tne  dispersion  relation  for  both  configurations.  The  two  solutions  ai’e  the  surface  or  brew star 
modes  and  the  virtual  modes.  The  characteristics  of  both  modes  are  discusseu. 

The  effect  oi  electronic  damping  upon  the  dispersion  curves  for  AG  which  exhibits  lov  elec- 
tronic damping  and  Au  which  exhibits  moderate  electronic  damping  is  demonstrated.  Finally  the 
perturbing  effect  of  the  dielectric  (referred  herein  as  dielectric  shift)  upon  the  displacement 
of  the  dispersion  curves  to  higher  wave  number  for  the  PAM  configuration  for  Ag  and  for  the  P31S 
configuration  for  Au  is  shown.  . 
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ABSTRACT 

The  dispersion  curves  for  surface  plasma  waves  (SPW)  in  Ag  have  been 
determined  from  calculated  reflectivity  minima  as  exhibited  by  attenuated  total 
reflection  (ATR  spectra)  for  the  prism-air-metal  (PAM)  configuration  and  from 
the  direct  calculation  of  the  dispersion  relation  for  the  same  configuration. 

The  dispersion  curves  for  Au  have  been  determined  by  measuring  the  ATR  spectra 
for  the  prism-dielectric-metal  (PDM)  configuration,  by  calculating  the  ATR 
spectra  from  published  optical  constants  and  from  the  direct  calculation  cf 
the  dispersion  relation  for  the  PDM  configuration.  We  have  found  two  general 
types  of  solutions  from  the  direct  calculation  of  the  dispersion  relation  for 
both  configurations.  The  two  solutions  are  the  surface  or  Brewster  modes  and  the 
virtual  modes.  The  characteristics  of  both  modes  are  discussed. 


The  effect  of  electronic  damping  upon  the  dispersion  curves  for  Ag 
which  exhibits  low  electronic  damping  and  Au  which  exhibits  moderate  electronic 
damping  is  demonstrated.  Finally  the  perturbing  effect  of  the  dielectric  (referred 
herein  as  dielectric  shift)  upon  the  displacement  of  the  dispersion  curves  to 
higher  wave  number  for  the  PAM  configuration  for  Ag  and  for  the  PDM  configuration 
for  Au  is  shown. 
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I - INTRODUCTION 
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Many  articles  [1-12]  have  been  published  concerning  the  determination 
of  the  surface  plasma  wave  (SPW)  dispersion  relation  for  metals,  both  theoretical 
and  experimental,  from  the  method  of  attenuated  total  reflection  (ATR)  [)].  The 
majority  of  attention  has  been  given  to  the  noble  metal  Ag  because  Ag  exhibits 
low  electronic  damping  and  therefore  large  dispersion.  Only  brief  attention  has 
been  given  to  metals  with  moderate  electronic  damping  such  as  Au  [9,10]  or  to 
metals  with  high  electronic  damping  such  as  nickel  [11,12]. 

In  this  study  the  SPW  dispersion  curves  for  Ag  as  determined  from 
calculated  ATR  spectra  for  the  PAM  configuration  and  from  the  direct  calculation 
of  the  dispersion  relation  for  the  same  configuration  are  presented.  The.  SPW 
dispersion  curves  for  Au  have  been  determined  from  experimental  ATR  spectra, 
from  calculated  ATR  spectra  and  from  the  direct  calculation  of  the  dispersion 
relation  for  the  PDM  configuration.  To  our  knowledge  this  is  the  first  reported 
direct  calculation  of  the  dispersion  relation  for  either  configuration.  We 
have  found  two  general  solutions  from  the  direct  calculation  of  the  dispersion 
relation  for  both  configurations.  They  are  the  surface  or  Brewster  modes  and  the 
virtual  modes  as  proposed  by  Kliewer  and  Fuchs  [13,1  A].  Also  we  have  found  that 
the  presence  of  the  dielectric  (prism)  in  proximity  to  the  air-metal  surface  for 
the  PAM  configuration  and  to  the  dielcctric-metal  surface  for  the  PDM  configuration 
causes  a shift  in  the  respective  dispersion  curves  to  higher  wave  numbers.  We 
refer  to  this  effect  as  the  "dielectric  shift".  In  effect  the  dielectric  shift 
is  a direct  result  of  the  inclusion  of  the  reflection  and  refraction  of  the  SPW 
by  the  second  boundary  at  a finite  distance  form  the  metal  surface. 

In  particular  we  give  attention  to  the  effect  of  the  dielectric  shift 
and  intrinsic  or  electronic  damping  upon  the  overall  shape  of  the  dispersion  curve 

as  determined  from  ATR  spectra  or  a direct  calculation  of  the  dispersion  relation 
as  derived  from  Maxwell's  equations.  Electronic  damping  is  directly 


3. 


proportional  to  the  magnitude  of  the  imaginary  part  of  the  complex  dielectric 
function  of  the  metal.  In  general,  the  larger  the  electronic  damping  the  smaller 
the  SPW  dispersion. 

The  original  configuration  of  the  ATR  method  for  inducing  and  detecting 
SPW  as  proposed  by  Otto  [1]  is  shown  in  Fig.  1.  An  alternate  arrangement  [2] 
in  which  the  air  film  is  replaced  by  a thin  film  of  the  active'  medium,  i.e.,  the 
metal,  is  shown  in  Fig.  2.  We  are  concerned  only  with  the  original  configuration 
in  this  study  ; the  second  arrangement  will  be  discussed  in  a subsequent  article. 

In  section  II  the  two-surface  dispersion  relation  is  presented  as 
derived  from  the  electromagnetic  wave  equation.  Comparison  is  made  with  published 
accounts.  A discussion  of  the  ATR  reflectivity  equation  is  also  given.  The  cal- 
culated dispersion  curves  from  the  two-surface  configuration  and  from  ATI!  spectra 
are  given  in  section  III  with  experimental  data.  Also  given  in  section  III  is 
an  explicit  evaluation  of  the  dielectric  shift.  A discussion  of  the  results 
obtained  is  presented  in  section  IV. 


II  - THEORY 

A)  Dispersion  relations 

The  SPW  is  ar.  electromagnetic  wave  propagating  as  a plane  wave  parallel 
to  the  interface  of  a metal  (active  medium)  and  a dielectric  but  decaying  expo- 
nentially perpendicular  to  the  interface.  The  presence  of  a second  dielectric  at 
a distance  t from  the  interface  modifies  the  propagation  behavior  of  the  SPW  [1). 
We  show  in  Fig.  3 the  electric  fields  in  the  three  media  configuration. 
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The  three  media  system  consists  of  a semi-infinite  dielectric  with 
dielectric  constant  cq,  a dielectric  layer  of  thickness  D with  dielectric 
constant  tj  (Cj  *=  I for  the  PAM  configuration)  and  a semi-infinite  medium  (m’tal 
or  semiconductor)  with  a complex  dielectric  constant  = c2r  + i c2i’  ^0=a^ 
limit  for  e2  is  assumed,  i.e., 

e2  (k,w)  -+  e2  (0,o) 

where  It  is  the  complex  wave  vector  and  u>  is  the  real  frequency  of  the  optical  field 

We  assume  that  an  evanescent  wave  has  excited  the  modes  of  collective 
oscillation  of  the  electron  density  at  the  surface  of  medium  2.  These  modes  cons- 
titute a SPW  propagating  parallel  to  the  surface  in  the  2 - direction  but  decaying 
exponentially  in  the  +x  direction.  The  method  of  "superposition  cf  waves"  is 
employed  to  account  for  waves  decaying  in  both  + x-directions  in  medium  1 . 

From  Maxwell's  equations,  the  wave  equation  that  must  be  satisfied  in 
all  three  media  is,  in  CGS  units, 

2 • 

(1)  V(V.E.)  - ”2E.  - ~ Ei  E.  - 0 i - 0,  !,  2 

c 

The  waves  in  the  three  media  are  transverse  magnetic  only.  SPVJ  in 
general  cannot  be  transverse  electric  waves.  In  the  three  medias  the  electric 
fields  are 
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(-  0,  i k)  E 


■^Op 


-anx 
c 0 


(2) 


1 


I.1  ■ <:  Sy-  o.  ‘ w i,p  *la'x 


h - <!:•  1 » h?  *v 


where  the  decay  constants  are 


2 ,2  u2 

a0  = k " k0  E0 


2 ,2  ,2 

“l  " k ‘ k0  E1 


2 2 2 
*2  - k - k0  C2 
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The  optical  wavelength  in  vacuum  is  X^.  The  time  and  z-dependcncy  factor  e1^2  ut‘ 
is  assumed  for  all  three  electric  fields,  Note  that  the  wave  number  ir.  the  z-direc— 
tion  is  simply  k. 


Applying  the  appropriate  boundary  conditions  at  the  two  interfaces  gives 
the  dispersion  relation  for  SPW  propagating  parallel  to  the  dielectric-metal  surfac. 
The  dispersion  relation  is 

al  e2  . (c0a1  + ela0)  + (e0al  • elV  e"2a,t  _ 

(->)  : — + — ■ ^2 — t " ® 

a2  El  (e0nl  + ela0)  " (E0al  " EIa0)  e~  • 


Note  that  the  second  term  in  Eq.  (5)  is  expressed  only  in  terms  of 
the  parameters  of  the  first  and  second  medium,  i.e.,  c^,  c^,  and  a(.  For 
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t very  large,  the  second  term  is  unity  and  the  dispersion  relation  for  a SPW 
propagating  parallel  to  a dielectric-metal  surface  obtains 


(6) 


2 C2r'l 

( £_L 

0 + e , 


When  t is  finite  the  deviation  of  the  second  term  in  Eq.  (5)  from 
unity  is  a direct  measure  of  the  influence  of  the  prism  upon  the  SPW  dispersion 
relation.  That  is,  it  is  possible  to  calculate  directly  the  effect  of  the  dielectri. 
shift  upon  the  SPW  dispersion  curve.  Eq.  (5)  was  programmed  for  computer  solution 
by  an  iteration  technique  to  yield  k = + i k. . The  solution  was  considered 

to  be  acceptable  when  the  magnitude  of  the  argument  of  the  left  side  of  Eq.  (5) 

-4 

was  less  than  10 


We  have  found  that  the  solutions  to  the  dispersion  relation,  Eq.  (5), 
for  the  two-surface  configuration  fall  into  tv?o  classes.  The  first  class  is 
characterized  by  the  decay  constant  cxq  in  the  pris’m  having  complex  values  with 
the  imaginary  part  positive.  The  real  part  of  is  generally  small  in  magnitude 
compared  to  the  imaginary  part  except  for  very  small  gap  thicknesses  at  the  smaller 
values  of  photon  energy.  The  sign  of  the  real  part  of  oiq  is  positive  for  this  class 
of  solution  except  for  intermediate  values  of  the  gap  thickness  and  certain  photon 
energies  where  it  is  negative.  The  other  decay  constants  a.  and  are  also  ccmplc:: 
with  positive  real  parts,  but  the  imaginary  parts  are  sometimes  positive  and  some- 
times negative  depending  on  the  photon  energy  and  the  gap  thickness.  The  wave  vec- 
tor k is  complex  with  a positive  real  part.  For  some  gap  thicknesses,  particularly 
at  the  lower  photon  energies,  the  imaginary  part  of  k is  found  to  be  negative, 
although  in  the  limits  of  infinite  gap  and  zero  gap,  it  is  always  positive.  When 
the  real  part  of  is  positive,  this  first  class  of  solution  corresponds  to  sur- 
face modes  of  the  type  discussed  by  Kliewer  and  Fuchs  [13,14]  in  which  the  fields 
decay  exponentially  into  the  inactive  medium. 
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Let  us  examine  the  components  of  the  real  part  of  the  Poynting  vector 
i>  in  the  prism  which  are  perpendicular  and  parallel,  respectively,  to  the  inter- 
faces : 


(7) 
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Bit 
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We  see  that  the  component  perpendicular  to  the  interfaces  is  proportional  to  a ^ , 
the  imaginary  part  of  Oq.  When  is  positive,  as  it  is  for  the  case  under  consi- 
deration, there  is  a flow  of  energy  in  the  prism  toward  the  active  medium.  In  this 
sense,  the  mode  is  a Brewster  mode  [20],  since  there  is  a single  electromagnetic 
wave  in  the  prism,  and  it  is  propagating  toward  the  active  medium. 

Let  us  now  turn  to  the  second  class  of  solutions  to  Eq.'  (5).  This  class 
is  characterized  by  the  decay  constant  having  real  and  imaginary  parts  which 
are  both  negative.  The  decay  constant  <Xj  has  positive  real  and  imaginary  parts 
for  the  cases  studied,  whereas  has  a positive  real  part,  but  either  a positive  o: 
negative  imaginary  part  depending  on  the  gap  thickness  and  photon  energy.  The  wave 
vector  k has  its  real  and  imaginary  parts  both  positive  for  all  cases  considered.  T! 
values  of  for  this  second  class  of  solution  corresponds  to  the  virtual  modes  of 

Kliewer  and  Fuchs  [13,14],  The  fields  increase  exponentially  going  away  from  the 
interfaces  into  the  prism,  and  energy  is  transported  away  from  the  interfaces. 

B)  Reflectivity  equations 

The  reflectivity  equations  required  to  determine  the  ATR  spectra  for 
the  prism-dielectric-mctal  system  as  shown  in  Fig.  1 have  been  discussed  exten- 


i 
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sively  in  the  literature  [1-12].  Therefore  only  the  final  equations  will  be 
presented  here. 


The  simplest  set  of  equations  for  calculating  the  ATR  spectra  is  given 
in  terms  of  a total  reflection  coefficient  R [15]  where 


(9) 


R = 


Rj  + R2  e 


-i  -2A 


1 + R!  R2  e 


-i  2A 


2 

The  reflectivity  is  simply  |r|  . 

Rj  is  the  Fresnel  reflection  coefficient  at  the  prism-dielectric  interface,  k2  is 
the  reflection  coefficient  at  the  dielectric-metal  interface  and  A is  the  phase 
factor  proportional  to  the  thickness  t of  the  air  gap.  The  Fresnel  reflection  coef- 
ficient Rj  must  be  modified  to  yield  total  internal  reflection  at  the  prism-dielec- 
tric surface  through  the  application  of  Snell's  law  [15]. 


An  alternate  form  of  the  total  reflection  coefficient  in  terns  of  the 
SPW  wave  vector  k and  the  decay  constants  an,  a.  and  a0  of  the  three  media  has 
been  given  by  Otto  [6,7].  A third  set  of  equations  has  been  presented  by  Kaplan 
et  al.  [16]  in  terms  of  scattering  matrices.  Both  Otto  and  Kaplan  et  al.  suggest 
that  setting  the  denominator  of  their  respective  total  reflection  coefficients  to 
zero  yields  the  SPW  dispersion  relation  for  a dielectric-dielectric-metal 
system.  Expressions  equivalent  to  Eq.  (5)  were  obtained  by  so  doing.  This  also 
served  as  a check  on  the  validity  of  the  model  as  shown  in  Fig.  3. 

An  equivalence  can  be  shown  between  Eq . (9)  and  the  total  reflection 
coefficient  as  presented  by  Kaplan  et  al.  [16]  by  using  the  decay  constants  ex- 
pressed in  terms  of  the  individual  bulk  dispersion  laws  for  each  medium.  The 
computer  calculations  of  Eq.  (9)  were  tested  by  programming  also  the  total 


reflection  coefficient  from  Kaplan  el  al.  [16].  Identical  results  were  obtained. 
All  calculations  of  Eq.  (9)  were  performed  by  varying  the  angle  of  incidence  by 
0.1  degree  increments  with  the  photon  energy  constant.  Published  optical  constants 
[17]  were  used  in  all  calculations  reported  herein.  The  real  part  of  the  surface 
plasma  wave  number  is  given  by  the  relation 

(10)  kr  «=  kQ  NQ  sin  0Q 

where  0^  is  the  angle  of  incidence  for  the  reflectivity  minimum  for  a constant 
photon  energy  and  Nq  is  the  index  of  refraction  of  the  prism. 


Ill  - RESULTS 

In  the  first  part  of  this  section  we  present  the  SPW  dispersion  curves 
for  Ag  and  Au  from  a direct  calculation  of  the  dispersion  relation  for  an  air- 
metal  surface,  i.e.,  Eq.  (C)  to  illustrate  the  dominating  effect  of  electronic 
damping.  In  the  second  part  the  SPW  dispersion  curves  for  Ag  for  the  PAM  confi- 
guration as  calculated  from  Eq.  (5)  and  Eq . (9)  are  presented.  Finally  the  SPW 
dispersion  curves  for  Au  for  the  FDM  configuration  as  calculated  from  Eqs.  (5) 
and  (9)  and  from  experimental  results  are  shown.  Experimental  results  from  lite- 
rature arc  also  included  where  possible. 

The  dielectric  constants  for  Ag  and  Au  are  plotted  as  a function  of 
photon  energy  in  Fig.  4.  The  dispersion  curves  for  Ag  and  Au  as  calculated  from 
Eq.  (6)  are  plotted  in  Fig.  5 for  c.  » 1.  It  is  instructive  to  compare  the  magni- 
tude and  behavior  of  Cj  ■ + i Ej^  for  the  two  metals  for  E ■ 2.0  - 4.0  eV 

with  regard  to  the  dispersion  relations.  For  Ag  J c 2^ I K<  lc2rl  ^or  E ^ 3*7  eV* 
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At  E = 3.6  eV  C2r  goes  through  - 1.0.  Then  the  denominator  on  the  right  side 
of  Eq.  (6)  becomes  very  small  and  thus  k ecoi.ies  large.  This  behavior  accounts 
for  the  large  dispersion' for  Ag  as  shown  in  Fig.  5.  The  dielectric  function  of 

Au  does  not  exhibit  this  type  of  behavior  and  therefore  Au  does  not  show  large 

« 

dispersion. 

SILVER 

The  virtual  mode  (a^  negative)  dispersion  curves  for  Ag  as  calculated 
from  Eq.  (5)  for  gap  thicknesses  t = 0.2  pm  and  0.075  pm  are  plotted  in  Fig.  6. 

The  latter  thickness  was  chosen  to  demonstrate  overcoupling.  The  solid  line  is  the 
dispersion  curve  for  t = ».  The  dielectric  (prism)  was  Al^O^.  The  dispersion  of 
the  index  of  refraction  for  AljO^  [18]  was  included  in  all  calculations.  The 
dispersion  curve  for  t = 0.2  pm  is  almost  identical  to  that  for  t = except 
in  the  photon  energy  range  of  2,0  to  3.0  eV  where  | | >>  1.  For  t «=  0,075  pm, 

the  dispersion  curve  exhibits  larger  wave  numbers  and  also  a larger  maximum  wave 
number . 


The  surface  mode  (a^  positive)  curves  for  Ag  as  calculated  from  F.q.  (5) 
for  gap  thicknesses  t «=  0.2  pm  and  0.075  pm  are  plotted  in  Fig.  7.  The  curves  are 
very  similar  to  those  for  the  virtual  modes  - in  fact,  for  t = 0.2  pm  they  are 
essentially  indistinguishable.  For  t «=  0.075  pm  the  surface  mode  curve  lies  slight! 
to  the  left  of  the  virtual  mode  curve  except  at  he  highest  photon  energies. 

The  direct  calculation  of  the  dielectric  shift  due  to  the  proximity  of 
the  prism,  i.e.,  the  second  term  on  the  left  hand  side  of  Eq.  (5),  for  t = 0.2  pm 
and  t **  0.075  pm  for  both  the  virtual  mode  and  the  surface  mode  are  shown  in 


Fig.  8.  We  see  that  the  deviation  of  the  dielectric  shift  from  unity  is  larger 
for  the  virtual  modes  than  for  the  surface  modes. 

The  dispersion  curves  for  Ag  as  determined  from  the  reflectivity  minima 
in  calculated  ATR  spectra  are  plotted  in  Fig.  9 for  air  gap  thicknesses  t = 0.2  pm 

and  0.075  pm.  The  dispersion  curve  for  t = «*  is  plotted  for  reference.  Two  dielec- 

tric light  lines  (0^  * 44.8°  and  60.5°)  are  also  shown.  For  t **  0.2  pm,  reflectivit 
minima  can  be  determined  for  photon  energies  3.5  eV.  For  t = 0.075  pm,  minima  can 
be  determined  up  to  3.6  eV  in  the  vicinity  of  which  backbending  occurs.  For 
t = 0.2  pm,  the  dispersion  curve  from  reflectivity  minima  is  practically  the  same 
as  those  for  both  the  virtual  modes  and  the  surface  modes.  For  t « 0.075  pm, 
however,  the  surface  mode  curve  seems  to  agree  slightly  better  with  the  reflecti- 
vity minima  curve  than  does  the  virtual  mode  curve. 

GOLD 

The  PDM  configuration  for  Au  was  fabricated  by  vacuum  depositing  on 
the  base  of  a 60°  heavy  dense  flint  prism  first  a thin  film  of  MgF2  (t  = 0.065  um) 

and  then  a thick  film  of  Au  (0.2  pm)  at  10  ^ torr.  The  ATR  spectra  were  measured 

in  the  photon  energy  range  of  1.0  - 2.6  eV.  Optical  coupling,  i.e.,  a reflectivity 
minimum,  could  -not  be  measured  above  E = 2.6  eV  because  of  the  inherent  angle  of 
incidence  limitation  in  the  ref lectometer . The  dispersion  curve  as  determined 
from  the  ATR  spectra  is  show:;  in  Fig.  10.  Also  the  dispersion  curve  from  the 
calculated  ATR  spectra  is  shown.  The  dispersion  of  the  indices  of  refraction  [19]  o 
the  two  dielectrics  was  included  in  all  calculations.  For  reference  the  dispersion 
curve  for  a dielectric  (MgF2)-metal  configuration  as  calculated  from  Eq.  (6)  is 
given.  Curve  fitting  of  £2  to  the  experimental  reflectivity  curves,  i.e.,  the 
ATI!  spectra  was  not  attempted  because  of 
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the  close  agreement  between  the  two  dispersion  curves. 


The  virtual  mode  (a^  negative)  dispersion  curve  and  the  surface  Mode 
( <oi  positive)  dispersion  curve  for  Au  for  the  PDM  configuration  (t  - 0.065  pm  for 
MgFj)  as  calculated  from  Eq.  (5)  are  presented  in  Fig.  11.  For  reference  the  dis- 
persion curve  from  the  calculated  ATR  spectra  is  also  shown.  The  surface  mode, 
dispersion  curve  is  a better  approximation  to  the  ATR  dispersion  curve  than  the 
virtual  mode  dispersion  curve  at  all  photon  energies.  Note  that  all  five  dispersio: 
curves  exhibit  the  same  general  shape,  particularly  in  the  photon  range  of  2.2  - 
2.6  eV.  This  illustrates  the  dominating  effect  of  the  moderate  electronic  damping 
of  Au. 
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The  real  and  imaginary  parts  of  the  SPW  number  for  Au  are  listed  in 
Tables  1 and  2 for  E = 1.14  eV  and  E • 2.26  eV  respectively  for  the  virtual  mode, 
the  surface  mode  and  from  ATR  spectra  for  several  dielectric  film  thicknesses. 

The  imaginary  part  was  calculated  from  the  ATR  spectra  using  the  relation 

(11)  k.  » k No  cos0„  A0 

1 O M 

where  0^  is  the  angle  of  incidence  at  the  reflectivity  minimum  and  A 0 is  the 
half-width  at  half -minimum  of  the  reflectivity  curve.  The  two  photon  energies 
were  selected  to  illustrate  the  difference  in  values  of  k close  to  the  dielectric 
, light  line,  i.e.,  at  E - 1.14  eV  and  far  from  the  dielectric  light  line,  i.o., 
at  E"“  2.26  eV.As  shown  in  Table  1 detectable  ATR  coupling  begins  at  t ■ 2.25  pm. 
At  t ■ 2 pm,  k for  both  the  surface  mode  and  virtual  mode  agree  with  k from  ATR 
spectra.  However  at  t ■ 1 ira  k.  for  the  surface  mode  is  negative  and  becomes  more 
negative  as  t decreases  to  t ■ 0.065  pm.  At  t ■*  0.065  pm  k^  is  smaller  in  magr.itud 
as  is  k^for  the  virtual  mode  and  ATR  spectra.  This  decrease  in  magnitude  indi- 
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cates  the  strong  overcoupling  that  occurs  between  t » 0.1  |ir  and  t = 0.065  urn. 

Detectable  ATE  coupling  begins  at  t = 0.25  um  for  E = 2.26  eV  as  shewn 
in  Table  2.  At  this  dielectric  film  thickness  is  larger  for  both  modes  than 
for  ATR  spectra  but  k^  is  smaller  for  both  modes  than  for  ATE  spectra.  For 
the  surface  mode  decreases  monotonically  but  dues  not  become  negative.  The 
surface  mode  is  a better  aproximation  to  the  ATR  spectra  than  the  virtual  mode  as 
can  be  seen  by  comparing  the  values  of  k^  for  all  film  thicknesses. 

The  direct  calculation  of  the  dielectric  shift  from  Eq . (5)  is  shown 
in  Fig.  12  as  a function  of  photon  energy  for  both  the  virtual  mode  and  the 
surface  mode.  Note  the  large  deviation  from  unity  in  the  lower  photon  energy 
range  where  | ^2r I >:>  ^8*  4)*  In  Fig.  13  the  dielectric  shift  is  plotted 

as  a function  of  dielectric  film  thickness  for  E = 1.14  eV  and  2.26  eV  for  both 
the  virtual  and  surface  modes.  Note  that  the  deviation  from  unity  is  more  nonlinea- 
and  begins  at  greater  film  thickness  at  E ■ 1.14  eV  than  at  E - 2.26  eV. 

Also  the  dielectric  shift  is  much  larger  for  the  lower  photon  energy.  In  both 
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Figs.  12  and  13  the  dielectric  shift  is  less  for  the  surface  mode  than  for 
the  virtual  mode. 


IV  - DISCUSSION 

We  have  found  that  there  are  two  general  types  of  solutions  to  the 
dispersion  equation  for  SPW,  Eq.  (5).  These  are  the  surface  or  Brewster  modes, 
in  which  the  decay  constant  in  the  prism  has  a positive  imaginary  part,  and 
the  virtual  modes,  in  which  has  a negative  imaginary  part.  Both  of  these  modes 
belong  to  a group  of  interface  electromagnetic  modes  which  involve  a single  elec- 
tromagnetic wave  in  the  prism  120].  The  surface  mode,  which  has  a positive  ima- 
ginary part  to  cxq,  involves  an  electromagnetic  wave  incident  on  the  prisrc-3ir 
interface,  but  no  reflected  wave.  This  mode  is  therefore  a Brewster-type  [20] 
mode.  The  virtual  mode,  on  the  other  hand,  has  a negative  imaginary  part  to 
and  therefore  involves  an  electromagnetic  wave  propagating  away  from  the  prism- 
air  interface,  but  no  incident  wave.  On  most  of  the  cases  studied  in  the  present 
work  the  imaginary  part  of  is  larger  in  magnitude  than  the  real  part  for  both 
the  surface  and  virtual  modes.  Consequently,  these  modes  possess  more  of  a propa- 
gating than  a decaying  character  in  the  prism. 

Let  us  discuss  for  a moment  the  virtual  mode  solutions.  The  real  part 
of  cXq  is  negative,  so  the  fields  increase  exponentially  going  away  from  the 
interfaces  into  the  prism,  and  energy  is  transported  away  from  the  interface.  As 
pointed  out  by  Kliever  and  Fuchs  [13]  one  should  introduce  a complex  frequency 
which  will  lead  to  a temporal  decay  of  the  wave  and  which  will  compensate  for 
the  exponential  growth  with  distance  into  the  prism.  Another  point  which  should 
be  mentioned  is  the  increase  in  value  of  ac  the  air  gap  thickness  is  decreased 
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from  large  values.  As  the  coupling  to  the  prism  increases  through  the  reduction 
in  gap  thickness,  the  rate  of  energy  transport  away  from  the  interfaces  increases, 
and  the  wave  must  therefore  decay  more  rapidly  in  the  direction  of  propagation. 
Although  this  behavior  was  observed  for  most  values  of  photon  energy,  it  was  not 
found  at  the  highest  photon  energies  where  backbending  occurs. 

Now  let  us  turn  to  the  surface  or  Brewster  modes.  In  this  case  the  value 
of  k^  decreases  as  the  air  gap  thickness  decreases  from  largv  values.  The  flow 
of  energy  in  the  prism  is  toward  the  interface  and  increases  as  the  gap  thickness 
decreases.  This  flow  of  energy  into  the  interface  region  tends  to  offset  the  decay 
of  the  mode  due  to  the  intrinsic  damping  in  the  active  medium  ; hence,  k^  decreases 
For  some  gap  thicknesses  and  some  photon  energies,  the  value  of  k^  is  negative.  On 
the  face  of  it,  this  would  imply  amplification  of  the  wave.  However,  there  is 
no  mechanism  available  for  amplification.  A possible  way  of  resolving  this 
situation  is  to  invoke  a complex  frequency  so  that  the  temporal  d.ecay  of  the  wave 
thereby  introduced  will  offset  the  spatial  growth  of  the  wave. 
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TABLES 


TABLE  1 - The  real  and  imaginary  parts  of  k (-  kf  + i k^)  for  Au  are  listed  for 

the  surface  mode,  virtual  mode  and  from  ATR  spectra  as  a function  of 

dielectric  film  thickness  t at  E ■ 1.14  eV. 

TABLE  2 - The  real  and  ima£ir.»r>  is  of  k (■  k ♦ i k^)  for  Au  are  listed  for 

the  surface  mode.  virtual  mode  and  from  ATR  spectra  as  a function  of 

dielectric  film  thickness  t at  E ■ 2.26  eV. 
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TABLE  ! 
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SURFACE  MODE 

VIRTUAL  MODE 
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ATR  SPECTRA  ) 

) 

( 

) 

( t 
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k. 

k 

k. 

k 

k.  ) 
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r 
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1 

r 

1 ) 

( 

) 

(2  pm 

0.808 

0.001 

0.808 

0.001 

0.808 

0.001  ) 

( 

) 

( 1 Pm 

0.806 

- 0.001 

0.806 

0.004 

0.806 

0.004  ) 

( 

• 

) 

(0.5  pm 

0.801 

- 0.021 

0.802 

0.025 

0.798 

0.023  ) 

( 

) 

0 

ro 

1 

0.840 

- 0.101 

0.843 

0.105 

0.816 

0.103  ) 

( 

) 

(0.1  pm 

0.959 

- 0.172 

0.964 

0.178 

0.898 

0.156  ) 

( 

) 

( 0.065pm 

1.106 

- 0.163 

1.116 

0.172 

0.997 

0.129  ) 

( 

1 

) 

( 

) 

ATR  coupling  begins  at  t = 2.25pm. 
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TABLE  2 
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SURFACE  MODE 

VIRTUAL  MODE 
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ATR  SPECTRA  ) 
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l 
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1 ) 

( 'A 

) 

( 0.25m 

1.855 

0.128 

1 .866 

0.137 

1.769 

0.155  ) 

( 

) 

( 0.2  pin 

1 .856 

0.118 

1 .882 

0.146 

1.814 

0.143  ) 

( 

% 

- ) 

( 0. 1 pm 

1.945 

0.050 

2.021 

0.226 

1.939 

0.216  ) 

( 

) 

( 0.065pm 

2.085 

0.040 

2.445 

0.281 

2.082 

0.221  ) 

( 

) 

( 

) 

ATR  coupling  begins  at  t * 0.25pm 
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! Fig.  1 - Prism-dielectric-metal  system  for  ATR.  At  0 >0  (0  * angle  of 

j “ o c c 

J , • 

j incidence  required  for  total  reflection  at  the  Cq  : Cj  surface) 

i 

an  evanescent  wave  created  at  the  prism  : dielectric  surface  couples 

| through  the  dielectric  film  to  excite  the  SPW  at  the  dielectric  : 

] metal  surface. 

I 

! 

4 

4 

| Fig.  2 - Prism-metal-air  (PMA)  system  for  ATR..  At  0 >0  an  evanescent  wave 

i — “ o c 

j created  at  the  metal  : air  surface  excites  the  SPW  at  the  same  surface. 

\ 

4 

j 

\ Fig.  3 - Schematic  diagram  showing  the  configuration  of  the  electric  field',  in 

the  three  media  for  the  derivation  of  the  SPW  dispersion  relation. 

Complex  dielectric  function  Ej  * e^r  + ^e2i  versus  photon  energy  E 
for  Ag  and  Au  [17). 

SPW  dispersion  curve  for  gold  .( ) and  silver  ( ) for  a metal-vacuum 

surface  as  calculated  from  published  optical  data  [17]  using  equation 
(6)  versus  photon  energy  E.  The  vacuum  light  line  is  k^  = 2-n/X^. 

Virtual  mode  dispersion  curves  for  Ag  for  the  PAM  configuration  for 

t “ 0.2  pm  ( ) and  t « 0.075  pm  (—  — ) versus  photon  energy  E.  The 

dispersion  curve  calculated  from  Eq.  (6)  ( ) and  the  vacuum  light 

line  k are  also  shown, 
o 
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Fig.  7 - Surface  mode  dispersion  curves  for  Ag  for  the  PAM  configuration  for 

t **  0.2  pin  ( ) and  t = 0.075  pm  (— ) versus  photon  energy  E.  The 

dispersion  curve  calculated  from  Eq . (6)  ( ) and  the  vacuum  light 

line  k are  also  shown, 
o 

Fig ■ 8 - Dielectric  shift  versus  photon  energy  E for  Ag  for  the  virtual  mode 

( ) and  surface  mode  ( ) at  t R 0.2  pm  and  0.075  pm  for  the  PAM 

configuration. 

Fig,  9 - SPW  dispersion  curves  for  Ag  for  the  PAM  configuration  from  the  calcu- 
lated ATR  spectra  for  t = 0.2  pm  ( ) and  t ■ 0.075  pm  (— — ) versus 

photon  energy  E.  The  experimental  data  (x..x)  from  Otto  [1]  and  two 
prism  light  lines  (0q  * 44.8°,  60.5°)  are  shown.  The  dispersion  curve 
calculated  from  Eq.  (6)  is  also  shown  for  reference. 


Fig.  10  - SPW  dispersion  curves  for  Au  for  the  PDM  configuration  from  the 

experimental  (— ) and  the  calculated  ( ) ATR  spectra  for  t * 0.065  pm 

versus  photon  energy  E.  The  calculated  dispersion  curve  for  t ■ " 

( ),  i.e.,  from  Eq.  (6)  (Cj  for  MgF^)  and  the  dielectric  light  line 

are  also  shown. 


Fig.  11  - Virtual  mode  (— ) and  surface  mode  ( ) dispersion  curves  for  Au  for 

the  PDM  configuration  from  Eq.  (5)  for  t “ 0.065  pm  versus  photon  energy 
E.  The  SPW  dispersion  curve  from  the  calculated  ATR  spectra  and  the 
dielectric  light  line  k^  are  also  shown. 


21. 

Pip-  1?  ” Dielectric  shift  versus  photon  energy  for  Au  for  the  virtual  mode 

(•—•—)  and  surface  mode  ( ) for  the  PDM  configuration  at  D = 0.065  pm. 

Fi E . 13  - Dielectric  shift  versus  dielectric  film  thickness  for  Au  for  the  virtual 

mode  (-—•-)  and  surface  mode  ( ) for  the  PDM  configuration  at  the  photc. 

energies  E * 1.14  eV  and  2.26  eV. 
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